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Abstract 

Measurements of inclusive production of the A, S~ and H*(1530) baryons in 
two-photon collisions with the L3 detector at LEP are presented. The inclusive 
differential cross sections for A and H~ are measured as a function of the baryon 
transverse momentum, pt, and pseudo-rapidity, rj. The mean number of A, S~ and 
S*(1530) baryons per hadronic two-photon event is determined in the kinematic 
range 0.4 GeV < pt < 2.5 GeV, \ri\ < 1.2. Overall agreement with the theoretical 
models and Monte Carlo predictions is observed. A search for inclusive production 
of the pentaquark ^"''(1540) in two-photon collisions through the decay 6^ — pKg 
is also presented. No evidence for production of this state is found. 



Submitted to Eur. Phys. J. C 



1 Introduction 



Two-photon collisions are the main source of hadron production in high-energy e+c^ interac- 
tions at LEP, via the process e'^e" — > e"'"e~7*7* — >• e'^e~hadrons, for which the cross section 
is many orders of magnitude larger than the e''"e~ annihilation cross section. The outgoing 
electron and positron carry almost the full beam energy and their transverse momenta are 
usually so small that they escape undetected along the beam pipe. At the LEP energies con- 
sidered here, the negative four-momentum squared of the photons, Q^, has an average value 
of ~ 0.2 GeV^ Therefore, the photons may be considered as "quasi-real" . In the Vector 
Dominance Model (VDM), each virtual photon can fluctuate into a vector meson, thus initi- 
ating a strong interaction process with characteristics similar to hadron-hadron interactions. 
This process dominates in the "soft" interaction region, where hadrons arc produced with a 
low transverse momentum, pt. Hadrons with high pt are instead mainly produced by the QED 
process 77 qq (direct process) or by QCD processes originating from the partonic content 
of the photon (resolved processes) . 

Fragmentation mechanisms can be investigated in two-photon reactions. These processes 
are described phenomenologically. In the Lund string model [1], hadron production proceeds 
through the creation of quark- ant iquark and diquark-antidiquark pairs during string fragmen- 
tation. Mesons and baryons are formed by colorless quark-antiquark and quark-diquark combi- 
nations, respectively. An extension of this model, the "simple popcorn mechanism" [2], includes 
the possibility of producing an additional meson between baryon-antibaryon pairs. The relative 
rate of occurrence of the baryon-meson-antibaryon configuration is governed by the so called 
"popcorn parameter." Many other parameters must be tuned to reproduce the measured hadron 
production rate, such as the strange-quark suppression factor, the diquark-to-quark production 
ratio or the spin-1 diquark suppression factor. 

In a statistical model approach [3], hadronisation is described with a reduced number of 
free parameters. Particles are produced in a purely statistical way from a massive colour- 
less gas which, for hadron production in two-photon collisions, is completely specified by two 
parameters: the energy density (or temperature) and a strange-quark suppression factor. In 
this model, the yield of different particles depends only on their masses, spins and quantum 
numbers. 

The L3 Collaboration has previously measured inclusive 7r°, Kg and A production in quasi- 
real two-photon collisions for a centre-of-mass energy of the two interacting photons, W^^, 
greater than 5 GeV [4,5]. In this Article this study is extended to the S^, S*(1530) and fl~ 
baryons^). The data sample consists of a total integrated luminosity of 610 pb^^, collected with 
the L3 detector [6] at e+e" centre-of-mass energies a/s — 189 — 209 GeV, with a luminosity 
weighted average value (y/s) — 198 GeV. Inclusive strange-baryon production has been exten- 
sively studied in e"'"e~ annihilation processes [7], whereas in two-photon collisions only inclusive 
A production has been previously measured at lower by the TOPAZ Collaboration [8]. 

Heavy-baryon detection techniques are also applied to search for pentaquark production in 
two- photon interactions, using the decay channel 9'^ — > pKg. Since the first observation of a 
narrow resonance near 1540 MeV in the K+n mass spectrum [9], interpreted as the pentaquark 
9'^, experimental evidence for and against this new state has been accumulated [10]. Second- 
generation experiments have so far not confirmed the existence of this state [11]. No search in 
two-photon reactions has been performed yet. 

^^If not stated otherwise, the symbols A, S~ and fl~ refer to both the A, H~ and f2~ as well as A, and 
fi^ baryons. All charge-conjugate final-states are analysed. 
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2 Monte Carlo simulation 



The process e^e~ c^c^hadrons is modelled with the PYTHIA [12] and PHOJET [13] event 
generators with twice the statistics of the data. In PYTHIA, each photon is classified as direct, 
VDM or resolved, leading to six classes of two-photon events. A smooth transition between 
these classes is obtained by introducing a pt parameter to specify the boundaries. The SaS- 
ID parametrization is used for the photon parton density function [14]. Since both incoming 
photons are assumed to be real in the original program, PYTHIA is modified to generate a 
photon fiux according to the equivalent photon approximation [15] with an upper Q"^ cut at the 
mass squared of the rho meson. 

PHOJET is a general purpose Monte Carlo which describes hadron-hadron, photon-hadron 
and photon-photon collisions. It relies on the dual parton model combined with the QCD- 
improved parton model [16]. The pt distribution of the soft partons is matched to the one 
predicted by QCD to ensure a continuous transition between hard and soft processes. The two- 
photon luminosity function is calculated in the formahsm of Reference 15. The leading-order 
GRV parametrisation is used for the photon parton density function [17]. 

In both programs, matrix elements and hard scattering processes are calculated at the 
leading order and higher-order terms are approximated by a parton shower in the leading-log 
approximation. The fragmentation is performed within the Lund string fragmentation scheme 
as implemented in JETSET [12], which is also used to simulate the hadronisation process. 
JETSET parameters are tuned by using hadronic Z decays [18]. In particular, the strangeness- 
suppression factor is set to 0.3, whereas the parameter governing the extra suppression of 
strange quarks in diquarks is fixed to 0.4 and the diquark-to-quark production ratio to 0.10. 
The popcorn parameter is set to 0.5 and a value otsimz) = 0.12 is used for the strong coupling- 
constant. 

The following Monte Carlo generators are used to simulate the background processes: 
KK2f [19] for the annihilation process c+c~ KORALZ [20] for e"'"e~ — >■ T"'"r~(7); 

KORALW [21] for e+e" W+W" and DIAG36 [22] for e+e" ^ e+e-r+r". The response 
of the L3 detector is simulated using the GEANT [23] and GHEISHA [24] programs. Time- 
dependent detector inefficiencies, as monitored during each data-taking period, are included in 
the simulations. All simulated events are passed through the same reconstruction program as 
the data. 

3 Two-photon event selection 

Two-photon interaction events are mainly collected by track triggers [25], with a track pt 
threshold of about 150 MeV, and the calorimetric energy trigger [26]. The selection of e"'"e~ — > 
e^e~hadrons events [27] is based on information from the central tracking detectors and the 
electromagnetic and hadronic calorimeters. It consists of: 

• A multiplicity cut. To select hadronic final states, at least six objects must be detected, 
where an object can be a track satisfying minimal quality requirements or an isolated 
cluster in the BGO electromagnetic calorimeter of energy greater than 100 MeV. 

• Energy cuts. The total energy deposited in the calorimeters must be less than 40% of ^/s 
to suppress events from the e"'"e~ — > qq(7) and e^e^ — >• r"'"r~(7) processes. In addition, 
the total energy in the electromagnetic calorimeter is required to be greater than 500 MeV 
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to suppress beam-gas and beam-wall interactions and less than 50 GeV to remove events 
from the annihilation process e+e~ —>■ qq(7). 

• An anti-tag condition. Events containing a cluster in the luminosity monitor with an 
electromagnetic shower shape, and energy greater than 30 GeV, are excluded from the 
analysis. The luminosity monitor covers the polar angular region 31 mrad < 6 < 62 mrad 
on both sides of the detector. In addition, events with an electron or positron scattered 
above 62 mrad are rejected by the cut on calorimetric energy. 

• A mass cut. The mass of all visible particles, Wms, must be greater than 5 GeV to 
exclude the resonance region. In this calculation, the pion mass is attributed to tracks 
while isolated electromagnetic clusters are treated as massless. 

About 3 million hadronic events are selected by these criteria with an overall efficiency of 
45% for a two-photon centre-of-mass energy, W^-y, greater than 5 GeV. The background level 
of this sample is less than 1% and is mainly due to the e+e~ qq(7) and e"'"e~ — >■ e^e~r+r~ 
processes. The backgrounds from beam-gas and beam-wall interactions are negligible. 

4 Strange baryon selection 

The A, S~ and Q~ baryons are identified through the decays A p^r^, Att^ and — >■ 

AK^. Due to their long lifetimes, the combinatorial background can be strongly suppressed 
by selecting secondary vertices which are clearly displaced from the interaction point. The 
S*(1530) decays strongly into S~7r"'" and no corresponding displaced vertex can be observed. 

The S~, Q~ and S*(1530) final states are reconstructed by combining A candidates with 
pion or kaon candidates. The latter are defined as tracks formed by at least 30 hits in the 
central tracker out of a maximum of 62 and a. pt > 100 MeV. The probability of the pion or 
kaon hypothesis, based on the dE/dx measurement in the tracking chamber, is required to be 
greater than 0.01. 

4.1 A selection 

The selection procedure is unchanged from our previous publication [5] . The A identification is 
optimized to achieve a high efficiency and a good background suppression by selecting secondary 
decay vertices satisfying the following conditions: 

• The distance d\ in the transverse plane^-* between the secondary vertex and the e"'"e~ 
interaction point must be greater than 3 mm. 

• The angle between the sum pt vector of the two tracks and the direction in the 
transverse plane between the e"'"e~ interaction point and the secondary vertex must be 
less than 100 mrad. 

The distributions of these variables are presented in Figure^ Good agreement between data 
and Monte Carlo is observed. The proton is identified as the track with the largest momentum, 
an assignment shown by Monte Carlo to be correct more than 99% of the time. The dE/dx 

^^The transverse plane is defined as the plane transverse to the beam direction. 
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measurements in the central tracking chamber must be consistent with this assignment, a 
probabihty greater than 0.01 for both the proton and the pion candidates being required. 

The distribution of the corresponding invariant mass of the pvr system, m(p7r), is shown 
in Figure |2t. A clear A peak over a smooth background is visible, compatible with the A 
baryon mass, rriA = 1115.68 ± 0.01 MeV [7]. The resolution of the signal, about 3 MeV, is 
well reproduced by Monte Carlo simulation. The m(p7r) mass spectrum for the pvr" and pvT"*" 
combinations are shown in Figure Eb and Efc, respectively. The numbers of A baryons for 
different pt and \t]\ bins are given in Tables Q and El They are determined by a fit in which the 
signal is modelled by a Gaussian and the background by a fourth-order Chebyshev polynomial. 
Consistent values of the fitted mass and width are found for the different pt and |?7| bins. 

4.2 S selection 

The S~ baryons are reconstructed by combining A candidates with pion candidates. For the 
identification of S~ and Q~ baryons, different criteria are used to select secondary vertices 
produced by A decays. The distance is required to be greater than 5 mm and the angle 
less than 200 mrad, as shown in Figure EK and Eb- The other cuts are left unchanged. The 
distribution of the resulting invariant mass m(p7r) is displayed in Figure Et and shows a clear A 
peak. The 76000 pvr combinations that lie in the mass interval 1.105 GeV < m(p7r) < 1.125 GeV 
are retained. To reduce the combinatorial background, the following criteria are applied: 

• the distance of closest approach (DCA) in the transverse plane of the S~ decay pion to 
the e"'"e~ interaction point, 6?,^, must be greater than 1 mm. 

• the distance in the transverse plane between the Avr vertex and the e"'"e~ interaction point, 
d^, is required to be greater than d^- 

• the angle a= between the pt vector of the Avr combination and the direction in the trans- 
verse plane between the e"'"e~ interaction point and the Avr vertex has to be less than 100 
mrad. 

Distributions of the difference d\ — d~ and of the angle a~ are displayed in Figure IH and 
exhibit a good agreement with the Monte Carlo predictions. 

The distribution of the mass of the Att system, m(A7r), is displayed in Figure El for the Avr" 
and Avr"*" combinations, respectively. Figure El shows the m(A7r) spectrum for the different pt 
bins listed in Tabled A clear S~ peak is visible over a smooth background, compatible with the 
measured S~ mass, m=- = 1321.31 ± 0.13 MeV [7]. The resolution of m(A7r), about 7 MeV, is 
well reproduced by Monte Carlo simulation. The number of H~ baryons in each pt and |?7| bin is 
evaluated by means of a fit to the m(A7r) spectrum in the interval 1.26 GeV < m(p7r) < 1.4 GeV. 
The signal is modelled with a Gaussian and the background by a fourth-order Chebyshev 
polynomial. The results are listed in Tables ^ and El 

4.3 S*(1530) selection 

The S*(1530) baryons are identified by reconstructing the decays S*(1530) —>■ 'E~7i~^ and 
S*(1530) — >■ S 7r~. The S~ candidates with a mass of ±15 MeV around the nominal S~ 
mass are combined with pion candidates. As the latter are produced at the e"'"e~ interaction 
point, their DCA must be less than 5 mm. To compare the S*(1530) production to other strange 
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baryons, the measurement is restricted to the kinematical region: 0.4 GeV < pt < 2.5 GeV, 
|?7| < 1.2. The distribution of the invariant mass of the Svr system, m(S7r), is shown in Figure 
for opposite-charge (S~7r"*" and 'E'^rc') and same-charge (H^tt^ and H^7r+) combinations. 
The number of same-charge combinations is normahzed to that of opposite-charge combina- 
tions in the region m(H7r) > 1.7 GeV. An excess corresponding to the H*(1530) is observed in 
the opposite-charge spectrum close to the nominal mass mH*(i53o) = 1531.80 ± 0.32 MeV [7]. 
The number of H*(1530) baryons is determined by means of a fit to the mass spectrum in the 
region 1.47 GeV < m(S7r) < 1.90 GeV, as shown in Figure Eb. The signal is modelled with a 
Gaussian and the background is parametrized by a threshold function of the form: 

a("^(S7r) - rnof exp c{m(^En) - mo) + (i(m(H^) - mo)' 
where a, b, c, d and mo are free parameters. The results are given in Table d 



4.4 Search for Q 

Since the topology of the — > AK^ decay is similar to that of Avr^, the selection 

criteria are identical except that kaon candidates are combined with A candidates instead of 
pion candidates. The corresponding m(AK) spectrum is displayed in Figure |H1 No signal is 
observed around the nominal mass of the Q~ baryon, tuq- = 1672.45 ± 0.29 MeV [7]. The 
expected signal, as predicted by PHOJET, is found to be almost negligible. 



5 Results and systematic uncertainties 

The cross sections for A, E~ and S*(1530) production are measured for W-y^ > 5 GeV, with 
a mean value (W^y) = 30 GeV, and a photon virtuality < 8 GeV^ with (Q^) ^ o.2 GeVl 
This kinematical region is defined by cuts at Monte Carlo generator level. 

The overall efficiencies for detecting A, S~ and S*(1530) baryons as a function of pt and 
\7]\ are listed in Tables ^ and |21 They include reconstruction and trigger efficiencies as well 
as branching fractions for the decays H*(1530) S^tt^, Att^ and A pvr^: 67%, 

100% and 64%, respectively [7]. The reconstruction efficiencies, which include effects of the 
acceptance and the selection cuts, are calculated with the PHOJET and PYTHIA Monte Carlo 
generators. As both generators reproduce well the shapes of the experimental distributions of 
hadronic two-photon production [27], the average selection efficiency is used. The track-trigger 
efficiency is calculated for each data-taking period by comparing the number of events accepted 
by the track triggers and the independant calorimetric-energy triggers. The efficiency of the 
higher level triggers is measured using prescaled events. The total trigger efficiency varies from 
82% for Pt < 0.4 GeV to 85% in the high pt region. 

The differential cross sections da/dpt, da/dpf and da/dlr^l for the reactions e"'"e~ e"'"e~AX 
and e"'"e~ — > e"'"e~S~X are given in Tableland Tabled respectively. The mean number of A, 
S~ and S*(1530) baryons per hadronic two-photon event is measured in the region 0.4 GeV < 
Pt < 2.5 GeV, \ri\ < 1.2 and W^^ > 5 GeV, as summarized in Table El Finally, the ratio of 
A to A and S~ to S baryons are determined from the respective mass spectra. They are 
found to be N(A)/N(A) = 0.99 ±0.04 and N(S^)/N(S-) = 0.96 it 0.14, where the uncertainties 
are statistical. These results are in agreement with the value of 1.0 expected for baryon pair 
production in two-photon reactions. 

The following sources of systematic uncertainties are considered: selection procedure, back- 
ground subtraction, limited Monte Carlo statistics, Monte Carlo modelling and the accuracy 
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of the trigger efficiency measurement. Their contributions to the cross section measurements 
are detailed in Table IH 

The uncertainties associated to the selection criteria are evaluated by changing the corre- 
sponding cuts and repeating the fitting procedure. The A selection uncertainty is dominated by 
the secondary vertex identification while the largest uncertainty for the H~ channel, about 9%, 
arises from the Avr vertex reconstruction. The uncertainty due to the e"'"e~ — >■ e^e~ hadrons 
event selection is less than 1%. 

The uncertainty due to background subtraction is assessed by using different background 
parameterizations and fit intervals in the fitting procedure. The Monte Carlo modelling uncer- 
tainty, taken as half the relative difference difference between PHOJET and PYTHIA, varies 
between 1% and 14% whereas the uncertainty associated to the limited Monte Carlo statistics 
varies from 1% for A to 18% for H*(1530). A systematic uncertainty of 2% is assigned to the 
determination of the trigger efficiency, which takes into account the determination procedure 
and time stability. 

6 Comparison with Monte Carlo and theoretical models 

The differential cross sections do"/dpj(e+e~ — >• e"'"e~AX) and do"/dpj(e+e~ e"'"e~S~X) for 
|?7| < 1.2 are presented in Figure The behaviour of the cross sections is well described 
by an exponential of the form y4exp(— ap^) in the region 0.16 GeV^ < pI < 1.69 GeV^ with 
a = 1.78 ± 0.16 GeV~^ and a = 1.68 ± 0.35 GeV~^ for A and S~ production, respectively. 
These two values are compatible with each other, as expected from the Lund string model. 
The region 1.69 GeV^ < < 6.25 GeV^ is better reproduced by a power law of the form Ap^^ 
with b = 5.2 ± 0.2 and b = 4.2 ± 0.9 for A and S~ baryons, respectively. 

The differential cross sections da/dpt{e'^e~ e"'"e~AX) and da/dpt{e'^e~ e"'"e~S~X) for 
It^I < 1.2 are displayed in FigurelHb- Phase-space suppression explains the lower values obtained 
in the ffist bins. The behaviour of the cross sections is well described by an exponential of the 
form Aexp{—pt/{pt)) in the region 0.75 GeV < pt < 2.5 GeV with a mean value {pt) = 368 ib 
17 MeV for A production and in the range 0.7 GeV < Pt < 2.5 GeV with {pt) = 472±98 MeV for 
S~ formation. These values are larger than those obtained for inclusive vr*^ and Kg production: 
(pt) ~ 230 MeV and {pt) ~ 290 MeV, respectively [4]. 

The data are compared to the PHOJET and PYTHIA Monte Carlo predictions in Figure 
ITUk . Both Monte Carlo programs fail to reproduce the shape and the normalization of the A 
cross section. A better agreement can be achieved by adjusting the width ag in the Gaussian 
Pt distribution for primary hadrons produced during the fragmentation. As shown in Figure 
ITUb . the predictions obtained with PYTHIA for A production using a value (Xg = 0.546 GeV 
reproduce much better the data than those with the default value ag = 0.411 GeV. On the 
other hand, an increase of the width of the pt distribution of quarks inside initial photons 
does not improve the description of the measured spectra. Similar trends are observed with 
PHOJET. 

The differential cross sections as a function of |?7| are displayed in Figure ITTk and ITTb . 
Both Monte Carlo programs describe well the almost-uniform \r]\ shape, while the size of the 
discrepancy on the absolute normalization depends on the pt range. 

The mean numbers of A, S~ and S*(1530) baryons per hadronic two-photon event are 
extrapolated from the measured phase-space to the full pt and |?7| range using the PYTHIA 
Monte Carlo with the adjusted value ag = 0.546 GeV. The uncertainty on the extrapolation 
factors propagates to the results as an additional systematic uncertainty. It arises mainly from 
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the modelling of the two-photon reaction and fragmentation processes. The contribution due to 
fragmentation, about 7%, is estimated using an independent fragmentation model [28] instead 
of the Lund string scheme to simulate baryon production as well as varying the parameters 
of these models. The two-photon modelling uncertainty, taken as the difference between the 
extrapolation factors estimated by PYTHIA and PHOJET, is found to be 3%. The results are 
summarized in Table El together with the predictions of PYTHIA and PHOJET. Both Monte 
Carlos overestimate the A and mean numbers, the discrepancy being less pronounced in 
PYTHIA than in PHOJET. A better agreement is obtained by decreasing the strange-quark 
suppression factor by a few percent. These extrapolated mean numbers can be converted to 
inclusive cross sections by multiplying them for the total two-photon hadronic cross section 
measured with the same detector as discussed in Reference [27]. 

The predictions of the thermodynamical model of Reference 3 are displayed in Figure IT^ as 
a function of the strangeness suppression factor 7^, using as energy density the value p = 0.4 [3]. 
Overall agreement is observed for a 7^ in the vicinity of 7^ = 0.7, which is similar to the value 
extracted from e"'"e~ annihilation events [29]. 

As suggested in Reference 30, the mean numbers (n) of octet and decuplet baryons can also 
be parametrized by a function of the form: 

{n) = a[(2J + l)/(2/ + 1)] exp{-bm^) 

where J, J, m are respectively the isospin, total angular momentum and mass of the baryons 
and a,b are free parameters. The ratios (2 J + 1)/(2J + 1) (n) are displayed in Figure IT^ 
together with the data measured in e"'"e~ coUisions at y/s = 10 GeV and i/i = 91 GeV [7]. A 
fit to our data gives b = 4.5 ± 0.9, a value compatible with the exponents b = 4.3 ± 0.6 and 
b = 4.0 ±0.1 obtained for e"*"e~ reactions at ^/s = 10 GeV and ^/s = 91 GeV, respectively. 
This provides evidence for the universality of fragmentation processes in two-photon and e"'"e~ 
reactions. 



7 Search for the 9^ pentaquark 

A 9^^^ Monte Carlo consisting of pentaquarks mixed with large samples of hadronic two- 
photon events is used to study the kinematics of the decay products and estimate the signal 
efficiency. In this Monte Carlo the pentaquarks are generated with a mass of 1.54 GeV and a 
width of 1 MeV. The transverse momentum and pseudo-rapidity distributions predicted by the 
PYTHIA event generator for inclusive S*(1530) production are taken to model the unknown 
distributions of the pentaquarks. As these distributions are very similar for different low mass 
baryons, the simulated acceptance is expected to be quite insensitive to the choice of a particular 
baryon. The — >• pKg decay is assumed to be isotropic. The Kgp and Kgp mass distributions, 
calculated at the generator level, show no evidence of a narrow peak close to 1540 MeV. This 
indicates that no reflections from known decay modes could generate a fake 9~^ signal. 

7.1 Event selection 

The 9~^ candidates are reconstructed using the decay 9^ pKg. Events are selected if they 
contain at least one track coming from the e"'"e~ interaction point and a Kg decaying into tt+tt" 
at a secondary vertex. The Kg selection proceeds as follows: 

■^^If not stated otherwise, the symbols 6*+ refers to 6*+ and 6 as both charge-conjugate final-states are 
analysed. 
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• Each track is required to have more than 12 hits out a maximum of 62 and pt > 100 MeV. 

• The dE/dx measurement of both pion candidates must be consistent with this hypothesis 
with a probabihty greater than 0.01. 

• The distance in the transverse plane between the secondary vertex and the e^e~ interac- 
tion point must be greater than 5 mm. 

• The angle between the sum pt vector of the two tracks and the direction in the transverse 
plane between the primary interaction point and the secondary vertex must be less than 
75 mrad. 

The distribution of the effective mass of the tttt system, m(7r7r), is shown in Figure IT^ A 
clear K3 peak is present over a smooth background. The spectrum is fitted with two Gaussian 
functions and a second order polynomial for the background. A central value of 497.6±0.1 MeV 
is obtained for the peak, which agrees with the expected value of 497.7 MeV [7]. To search for 
, the Kg candidates are selected within an interval of ±20 MeV around the central value of 
the peak, corresponding to 140000 Kg candidates with a purity of 69%. 

The proton and antiproton candidates are tracks with at least 30 hits out of a maximum 
of 62 and a DCA less than 5 mm. In addition, the dE/dx measurement in the central tracker 
must be consistent with the proton hypothesis with a probability greater than 0.05. To reject 
kaons and pions, the probability of each of these hypothesis must be less than 0.01. A total of 
7131 selected protons and 5340 antiprotons are selected. Their dE/dx measurement is shown 
in Figure El together with theoretical the expectations based on the Bethe-Bloch formula. The 
purity of the selected sample is greater than 96%. The excess of protons, due to secondary 
interactions between particles and the detector, is well reproduced by the Monte Carlo. 

The resulting pKg and pKg mass spectra, m(pKg), are shown in Figure IT^ and 116b respec- 
tively. Since no distinct structure due to the excess of protons is observed, the two spectra are 
combined in the following. 

7.2 Results 

An upper limit for the number of e"'"e~ e^e^0+X events at 95% confidence level is derived 
with a fit to the m(pKg) mass spectrum. To minimize the effect of the binning and the 
choice of the fitting interval, an unbinned maximum likelihood fit is performed in the range 
1.45 GeV < W^^ < 1.8 GeV. The likelihood function is written as: 

C = p- c/(m(pKO)) + (1 - p) ■ 6(m(pK0)) 

where p denotes the fraction of signal events inside the fitted region. The background, 6(m(pK0)), 
is parametrized by a fourth-order polynomial and the signal by a Gaussian distribution, g{m(^p^o-j) 
The resolution of the expected signal is determined from Monte Carlo to be 14.0 ± 0.6 MeV. 
The resulting fit is displayed in Figure IT7I and yields a fraction of signal events compatible with 
zero: p = —0.006 ± 0.007. Restricting the measurement to the physical region p > 0, an upper 
limit of 59.3 events at 95% confidence level is derived. 

The selection efficiency is estimated to be 1.0%, including a branching fraction Kg 
TT+TT^ of 69% [7]. The branching ratio 6^ pKg is set to 1/4, accounting for competition with 
both the K'*"n and KlP channels. The overall trigger efficiency is measured to be 83%. 
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The mean number of 9^ pentaquark per hadronic two-photon event extrapolated to the full 
phase-space for VF-y^ > 5 GeV is found to be less than 4.0 x 10~^ at 95% confidence level. The 
systematic uncertainty arising from the selection procedure is evaluated by varying the cuts 
and repeating the fitting procedure. The contribution of the Kg reconstruction is estimated to 
be 9% whereas an uncertainty of 4% is associated with proton identification. The uncertainty 
due to 9'^ modelling, determined by using the pt and pseudo-rapidity distributions of different 
baryons in the 9^ generation, is found to be 9%. The effect of a non-isotropic 9'^ — > pKg decay 
on the efficiency is estimated using an angular distribution of the form 1 — cos^ a, where a is 
the angle of the proton in the 9^ ccntrc-of-mass system. The associated systematic uncertainty 
is found to be 6%. Limited Monte Carlo statistics introduces an additional uncertainty of 
5%. The sum in quadrature of these contributions yields a total systematic uncertainty of 16%. 
Including this uncertainty in the determination of the number of 9'^ pentaquark per two-photon 
hadronic event gives 

{ne+) < 4.7 x lO'^ 

at 95% confidence level. This result is about four times greater than the mean S*(1530) 
multiplicity: a baryon with almost the same mass. This difference is mainly due to the stringent 
cuts apphcd to proton selection, resulting in a smaller selection efficiency than that of the 
S*(1530) baryon, as well as the low branching ratio ^+ — > pKg, considered as half of that of 
the decay S*(1530) — > A7r~7r+ — > p7r^7r~7r''". 

Using the two-photon cross-section for W^^j > 5 GcV, cr^^ = 397 nb [27], a 95% C.L. upper 
limit on the cross section 77 9'^X of 1.8 nb is obtained. This result is comparable to the 
upper limit obtained by the photoproduction experiment in the reaction 7p — > K^9^ [11]. 

As a cross check, several subsets of the selected pKg samples were investigated using tighter 
selection criteria. The Kg selection interval was reduced to ±10 MeV, the DCA of proton and 
antiprotons candidates was required to be less than 3 mm or secondary vertices compatible 
with the A mass hypothesis were rejected. Since strangeness is a conserved quantity in strong 
interactions, the presence of an additional kaon was also required. No significant 9~^ signal was 
observed in any of these samples. 

8 Conclusion 

The production of A, S~ and S*(1530) baryons in two-photon collisions is studied in the range 
0.4 GeV < Pt < 2.5 GcV, \r]\ < 1.2 and W^^ > 5 GcV. The shape of the differential cross 
section for A and production is relatively well reproduced by the PYTHIA and PHOJET 
Monte Carlo programs using parameters tuned at y/s = rriz although a better agreement 
can be obtained for A by increasing the width in the gaussian pt distribution for primary 
hadrons. The mean numbers of A, S~ and S*(1530) baryons per hadronic two-photon event 
are found to be slightly below the PYTHIA and PHOJET predictions but overall agreement 
with the thcrmodynamical model is observed. The comparison between measurements obtained 
in two-photon events and e"'"e~ annihilation processes provides evidence for the universality of 
fragmentation functions in both reactions. 

Finally, a search for the pentaquark ^"'"(1540) through the decay 9^ — > pKg for W^^ > 5 GeV 
is presented. No evidence for 9~^ production is found and a 95% confidence level upper limit 
on the mean number of 9^ per two-photon hadronic event at a level of four times the observed 
rate for the S*(1530) baryon is derived. 
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Pt 

(GeV) 


(Pt) 
(GeV) 


( GeV') 


Number of 
baryons 


Efficiency 
(%) 


da/ dpt 
(pb/ GeV) 


da/dpf 
(pb/ GeV') 


A baryon 


0.4- 


-0.6 


0.50 


0.25 


3412 ± 71 


10.2 ±0.1 


273.1 ± 5.7 ± 35.7 


273.1 ± 5.7 ± 35.7 


0.6- 


-0.8 


0.69 


0.48 


4408 ± 85 


13.7±0.1 


264.1 ± 5.1 ± 29.1 


188.7 ± 3.6 ± 20.8 


0.8- 


-1.0 


0.89 


0.79 


3420 ±81 


15.3 ±0.2 


183.3 ± 4.4 ± 15.3 


101.9 ± 2.4 ± 


8.5 


1.0- 


-1.3 


1.12 


1.27 


3201 ±87 


16.8 ±0.2 


103.9 ± 2.8 ± 


7.8 


45.2 ± 1.2 ± 


3.4 


1.3- 


-1.6 


1.43 


2.05 


1222 ±55 


17.7 ±0.4 


37.8 ± 1.7 ± 


3.1 


13.0 ± 0.6 ± 


1.1 


1.6- 


-2.0 


1.77 


3.14 


578 ± 41 


15.9 ±0.6 


15.0 ± 1.1 ± 


1.4 


4.2 ± 0.3 ± 


0.4 


2.0- 


-2.5 


2.21 


4.92 


292 ± 25 


17.2 ± 1.3 


5.6 ± 0.5 ± 


1.2 


1.2 ± 0.1 ± 


0.3 


S baryon 


0.4- 


-0.7 


0.55 


0.31 


70± 10 


3.4 ±0.2 


11.3 ± 1.7 ± 


1.7 


10.3 ± 1.5 ± 


1.6 


0.7- 


-1.0 


0.83 


0.70 


113 ± 12 


7.1 ±0.3 


8.7 ± 1.2 ± 


1.2 


5.1 ± 0.6 ± 


0.7 


1.0- 


-1.3 


1.13 


1.27 


83± 12 


9.7 ±0.8 


4.7 ± 0.8 ± 


0.8 


2.0 ± 0.3 ± 


0.3 


1.3- 


-2.5 


1.67 


2.88 


94± 19 


8.8 ±0.9 


1.5 ± 0.3 ± 


0.3 


0.4 ± 0.1 ± 


0.1 












*(1530) baryon 








0.4- 


-2.5 


0.82 


0.84 


56 ±13 


3.8 ±0.6 


1.4 ± 0.3 ± 


0.5 


0.5 ± 0.1 ± 


0.2 



Table 1: The average transverse momentum (pt), the average transverse momentum squared 
(pj), the number of baryons estimated by the fits to the A, and S*(1530) mass spectra, the 
detection efficiency and the differential cross sections da/dpt and dcr/dpf for inclusive A, E~ 
and S*(1530) production as a function of pt for \r]\ < 1.2. The first uncertainty on the cross 
sections is statistical and the second systematic. 



1^1 


(1^1) 


Number of 


Efficiency 


dcr/d ?7 








baryons 


(%) 


(pb) 


A baryon 0.4 GeV < Pt < 1.0 GeV 


0.0- 


-0.3 


0.15 


2953 ±72 


13.8 ±0.2 


58.6 ± 1.4 ± 3.4 


0.3- 


-0.6 


0.45 


2742 ± 63 


13.4 ±0.2 


56.1 ± 1.3 ± 3.2 


0.6- 


-0.9 


0.75 


2904 ± 70 


13.0 ±0.2 


61.0 ± 1.5 ± 4.1 


0.9- 


-1.2 


1.05 


2774 ± 89 


11.1±0.1 


68.0 ± 2.2 ± 8.4 


A baryon 1.0 GeV < Pt < 2.5 GeV 


0.0- 


-0.3 


0.15 


1458 ± 60 


18.8 ±0.5 


21.2 ± 0.9 ± 1.6 


0.3- 


-0.6 


0.45 


1411 ±62 


18.2 ±0.4 


21.2 ± 0.9 ± 1.6 


0.6- 


-0.9 


0.75 


1480 ± 62 


19.3 ±0.5 


20.9 ± 0.9 ± 1.6 


0.9- 


-1.2 


1.05 


1007 ±62 


14.3 ±0.4 


19.3 ± 1.2 ± 2.7 






=] baryon 0.4 


[ GeV <pt< 2.5 GeV 


0.0- 


-0.4 


0.20 


110 ± 14 


6.7 ±0.3 


6.9 ± 0.9 ± 1.0 


0.4- 


-0.8 


0.60 


114 ± 13 


6.1 ±0.2 


8.0 ± 1.0 ± 1.0 


0.8- 


-1.2 


1.01 


109 ± 15 


5.5 ±0.2 


7.6 ± 1.1 ± 1.2 



Table 2: The average pseudo-rapidity the number of baryons estimated by the fits to the 
A and S~ mass spectra, the detection efficiency and the differential cross sections dcr/dlr^l for 
inclusive A and S~ production as a function of |?7|. The first uncertainty on the cross sections 
is statistical and the second systematic. 
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Data 


PYTHIA 


PHOJET 




0.4 GeV <pt< 2.5 GeV, rj < 1.2 and W^^ > 5 GeV 


/r, \ 

\n=.(i530); 


(^i.o ± U.iJ X iU 

[ 1 .0 ± i.uj X iu 
(2.3 ± 1.0) X 10-^ 


[lAo ± U.Ui j X iU 

(Q fil J- n 1 1 ^ \x 1 n— 4 
(^o.Oi ± U.iij X iU 

(1.37 ±0.04) X 10-^ 


^^ on _i_ n m N ^ l n— 2 
(^i.oU ± U.Ui j X iU 

(^ii.O ± U.iUj X iU 

(1.91 ±0.03) X 10-^ 












Data 


PYTHIA 


PHOJET 




W^^ > 5 GeV 


{nA) 

(?T'S*(1530)) 


( 7.6 ±0.8) X 10-2 
( 3.7 ±0.5) X 10-3 
(11.8 ±5.2) X 10-^ 


(8.51 ±0.01) X 10-2 
(5.17 ±0.02) X 10-3 
(8.65 ± 0.09) X 10-^ 


(1.11 ±0.01) X 10-1 
(7.29 ± 0.09) X 10-3 
(12.4 ±0.10) X 10-3 



Table 3: The mean number of A, S" and S*(1530) baryons per hadronic two-photon event for 
the measured pt and I77I phase-space (top) and its extrapolation (bottom) to the full pt and 
I77I range. The uncertainty on data includes both the statistical and systematic contributions 
whereas the uncertainty on the Monte Carlo predictions is statistical. 



Pt 


Selection 


Background 


Monte Carlo 


Monte Carlo 


Total 


(GeV) 


criteria 


subtraction 


statistics 


modelling 




A baryon 


0.4-0.6 


4.2 


12.1 


0.9 


1.4 


13.1 


0.6-0.8 


4.2 


9.8 


0.9 


1.6 


11.0 


0.8-1.0 


4.2 


6.6 


1.2 


1.9 


8.4 


1.0-1.3 


4.2 


5.3 


1.4 


2.1 


7.5 


1.3-1.6 


4.2 


5.8 


2.4 


2.4 


8.2 


1.6-2.0 


4.2 


6.4 


3.5 


2.8 


9.1 


2.0-2.5 


4.2 


19.4 


7.5 


4.8 


21.8 


S baryon 


0.4-0.7 


10.3 


8.7 


5.8 


3.1 


15.1 


0.7-1.0 


10.3 


7.5 


4.8 


2.4 


14.0 


1.0-1.3 


10.3 


9.1 


8.3 


3.3 


16.5 


1.3-2.5 


10.3 


12.3 


10.4 


6.3 


19.5 


S*(1530) baryon 


0.4-2.5 


19.2 


24.6 


16.5 


13.6 


37.9 



Table 4: Systematic uncertainties in percent on the cross section of the e+e- e+e^AX, 
e"'"e- e+e-5-X and e+e- —>■ e"'"e-S*(1530)X processes due to selection criteria, background 
subtraction, hmited Monte Carlo statistics, Monte Carlo modeUing and trigger efficiency as a 
function of pt for I77I < 1.2. The total systematic uncertainty, taken as the quadratic sum of 
the different contributions, includes a constant contribution of 2% due to the determination of 
the trigger efficiency. 
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Figure 1: Distribution of the variables used to select secondary vertices from inclusive A pro- 
duction: a) the distance in the transverse plane between the secondary vertex and the e'''e~ 
interaction point, dA, and b) the angle between the sum pt vector of the two tracks and the di- 
rection in the transverse plane between the interaction point and the secondary vertex, a\. All 
other secondary vertex selection criteria are applied. The predictions of the PHOJET Monte 
Carlo are shown as the sohd lines and the contribution due to A baryons as the dashed lines. 
The Monte Carlo distributions are normahzed to the data luminosity. 
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Figure 2: The mass of the pvr system for the inclusive A selection for a) A and A candidates, b) A candidates and c) A candidates. The 
signal is modelled with two Gaussians and the background by a fourth-order Chebyshev polynomial. 
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Figure 3: Distribution of the variables used to select secondary vertices from A produced in E~ and fl~ decays: a) the distance in the 
transverse plane between the secondary vertex and the c+c^ interaction point, dA, and b) the angle between the sum pt vector of the 
two tracks and the direction in the transverse plane between the interaction point and the secondary vertex, a^. All other secondary 
vertex selection criteria are applied. The predictions of the PHOJET Monte Carlo are shown as the full line and the contribution due to 
A baryons as the dashed hne. c) The mass of the pn system for A baryons produced in S~ and Q~ decays. The signal is modelled with 
two Gaussians and the background by a fourth-order Chebyshev polynomial. Only the combinations in the ±10 MeV mass window 
indicated by the arrows are retained to reconstruct E~ and fl~ baryons. The Monte Carlo distributions are normalized to the data 
luminosity. 




Figure 4: Distribution of the variables used for the selection of S~ baryons: a) the difference 
dA — ds of the A and S~ vertex distances from the interaction point and b) the angle as 
between the pt vector of the Avr combination and the direction in the transverse plane between 
the primary interaction point and the Avr vertex. In each plot, all other selection criteria 
are applied. The predictions of the PHOJET Monte Carlo are shown as the full line and 
the contribution due to S~ baryons as the dashed line. The Monte Carlo distributions are 
normalized to the data luminosity. 
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Figure 5: The mass spectrum of the Att system for a) S~ and b) S candidates for 0.4 GeV < 
Pt < 2.5 GeV and \ri\ < 1.2. The signal is modelled with a Gaussian and the background by a 
fourth order Chebyshev polynomial. 
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Figure 6: The mass of the Att system for a) 0.4 GeV < pt < 0.7 GeV, b) 0.7 GeV < pt < 
1.0 GeV, c) 1.0 GeV < pt < 1.3 GeV and d) 1.3 GeV < pt < 2.5 GeV. The signal is modelled 
with a Gaussian and the background by a fourth-order Chebyshev polynomial. 
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Figure 7: a) The mass of the Svr system for opposite charge ('E~7r~^ and S tt^) and same-charge 
(S~7r~ and S^tt"'') combinations. The number of same-charge combinations is normahzed 
to that of opposite-charge combinations in the region msyr > 1-7 GeV. A signal consistent 
with S*(1530) production is observed, b) Fit to the mass spectrum of the opposite-charge 
combinations. The signal is modelled with a Gaussian and the background by a threshold 
function. 
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Figure 8: The mass of the AK system. No Q~ signal is observed around the mass of 1.67 GeV. 
The predictions of the PHOJET Monte Carlo are shown as the full line and the contribution 
due to Q~ baryons as the dashed histogram. The Monte Carlo distributions are normalized to 
the data luminosity. 
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Figure 9: a) Differential cross sections da/dpf for the e+e" e'^e^AX and e+c^ c+e~S~X 
processes for |?7| < 1.2. b) Differential cross sections da/dpt for the e+e" — > e+e^AX and 
e+e~ — > e"'"e~S~X processes for \r]\ < 1.2. The various fits are described in the text. 
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Figure 10: The differential cross section as a function of pt for A and S~ production for 

\r]\ < 1.2 compared with a) the predictions of the PHOJET and PYTHIA Monte Carlo, and 
b) the predictions of the PYTHIA Monte Carlo using the default value for the width (Tq in 
the Gaussian pt distribution on primary hadrons Ug = 0.411 GeV as well as the adjusted value 
ag = 0.55 GeV. 
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Figure 11: The differential cross section as a function of \ri\ for a) inclusive A production for 
0.4 GeV < pt < 1.0 GeV and 1.0 GeV < pt < 2.5 GeV and b) for inclusive S~ production for 
0.4 GeV < Pt < 2.5 GeV. The predictions of the PHOJET and PYTHIA programs are also 
shown. 
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Figure 12: The extrapolated mean multiplicities (n) of A, S~ and S*(1530) baryons per two- 
photon event for W^^ > 5 GeV (horizontal bars) compared to the predictions of the thermody- 
namical model as a function of the strangeness suppression factor, 7^, using as energy density 
the value p = 0.4 (symbols). 
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Figure 13: Weighted mean multiplicities (2 J + I)/ {21 + 1) (n) as a function of the square 
of the baryon masses measured in two-photon reactions (solid circles) together with the data 
obtained in e"'"e~ collisions at ^/s = 10 GeV (open squares) and ^/s = 91 GeV (open circles). 
The hnes correspond to the fits described in the text. 
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Figure 14: The effective mass of the tttt system. The signal is modelled with two Gaussian 
functions and the background by a second-order polynomial. About 140000 Kg candidates are 
found in a ±20 MeV window around the central value of the peak and used in the search for 
the 9'^ pentaquark. 
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Figure 15: The d£'/dx measurement as a function of the momentum for the selected a) protons 
and b) antiprotons, together with the calculations of the Bethe-Bloch formula for protons, kaons 
and pions. 




Figure 16: The mass of the pKg system, m(pKg), for a) protons and b) antiprotons together 
with the predictions of the PYTHIA Monte Carlo. The arrow indicates the position of the 
expected 6'^ signal, around 1.54 GeV. 
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Figure 17: The mass of the pKg system, m(pKg), for protons and antiprotons combined together 
with the fit result. The arrow indicates the position of the expected 0^ signal, around 1.54 GeV. 
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